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Abstract 

We have investigated a cosmological model with variable speed of light (c) , gravitational constant 
(G) and cosmological constant (A). The model is shown to solve the horizon, flatness and monopole 
problems of the early universe. We have found that with a certain variation of these parameters the 
model predicts a cosmic acceleration. The model also predicts that for a flat universe A vanishes in 
both radiation and matter phases. If the gravitational constant is allowed to increase then we might 
not need the existence of dark matter. 
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1 Introduction 

The idea of change fundamental constant of physics was first suggested by Dirac (Dirac, 1937, 1938). He 

postulated that the gravitational constant (G) decreases with time (t) as, G oc \ . It was thought such 

a variation would help understand the existence of very large numbers appearing when one compares 

atomic physics with cosmology. But very recently Moffat (1993), Albrecht and Magueijo (1999), Barrow 

(1999, 2003), Avelino & Martins (1999) have conjectured that if the speed of light falls at certain rate 

then the horizon, flatness and monopole problems of the standard model can be solved without recourse 
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to inflationary paradigm. However, the variation of the speed of light will have important consequences 
especially during the nucleosynthesis era. Besides this, the variation of the speed of light will require a 
critical revision of the theory of relativity which relies upon the constancy of speed of light. The implica- 
tion of this variation will be immense in the arena of thermodynamics, astrophysics and cosmology. This 
variation will also alter our basic principle and laws in physics, e.g., Lorentz invariance and general co- 
variance. It seems plausible, that the force of gravity can be integrated into the theory of electrodynamics 
if the velocity of light changes according to some cosmological law. 

We will postulate here in this letter that such a change is countered by a change in the cosmological (or 
gravitational) constant, in such a way that the usual energy conservation still holds. With this minimal 
change of cosmology in mind we investigate the consequences of this variation. 

2 The model 

Solution of Einstein field equations for a universe with cosmological constant A yields the two equations 

/ R\ 2 8itG P kc 2 Ac 2 



+ ^r, (i) 



and 



\RJ 3 R 2 3 

( R\ AttG, 3p. Ac 2 

The cosmic fluid is characterized by the equation of state 

p = ( 7 -l)pc 2 1<7<2. (3) 

It is assumed by Barrow (2003) that the Friedmann cosmology in the presence of time dependent speed 
of light c{t) is described by the same Friedmann equation as those of constant c. Therefore, in what 
follows we will consider c and A as functions of time (or scale factor). We employ the following ansatzs 
for A and c as 

c = c R n , A = A R m , (4) 
where n,m, cq, Aq are undetermined constants. With this prescription, eqs.(l) and (2) yield 

(* V = tfn-* + ^Aocg +2n 37 (5) 

\R 2n + 37 - 2 m + 2n + 3j ' w 
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where A is an integration constant. Now consider the following cases: 
If A ^ 0. We see from eq.(5) that whenever the term 

n<^(2-37) and m < -2 , (6) 

first term in LHS of eq.(5) becomes negligibly small for large R in comparison with the other two terms. 
If this happens the curvature term will be insignificantly small and the flatness problem is automatically 
solved, as well as the horizon and monopole problems. We further admit that if the universe is dominated 
by strings (7 = |), then there will be no flatness or horizon problems, as evident from eq.(5). Inflationary 
solution is known to require the condition 7 < — | implying an exotic equation of state to solve these 
problems. We remark here that our eqs.(5) and (6) generalize Barrow (2003) equations to include a 
variable A term. We arrive at eqs.(5) and (6) using Friedmann equations directly, in comparison with 
Barrow (2003) method which employs matter conservation law [see eq.(3)]. 
In the radiation dominated phase 7 = | so that eq.(6) reduces to 

n < -1 (7) 

which implies a decaying speed of light. 

3 Early acceleration 

If the universe was dominated by strings and that 0, then irrespective of the value of k (i.e., 0, 1, —1), 
one would have 



yRJ m + 2n + 37 

if A = 0. We see that the curvature term vanishes identically. It is evident from the above equation that, 
when Ao > one has 

m + 2n + 2 > (9) 
With this constraint eq.(8) can be easily solved to give 

R = Bt~^^ (10) 

where B= constant. The constraint in eq.(9) implies that 

R>0. (11) 



Thus one gets a power law inflation in comparison with the problematic inflationary solution which 
normally requires an existence of an unseen scalar particles. 

Now if A ^ but k = 0, A = 0, one has for radiation dominated (RD) phase the equation 

(l) 2 ^-, (12) 

which can be solved to give 

R = Dtf (13) 



where D = const. This is the usual Einstein-de Sitter solution for RD phase. 

4 Matter dominated (MD) phase 

For matter dominated phase one has 7 = 1. For a flat universe (k = 0) with A = 0, A ^ 0, one has 

(fY = ^-3 (14, 



which is solved to give 
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R = D'ts , (15) 



which is the usual Einstein-de Sitter solution. We observe that this cosmology is determined by the 
constants k, c, A and does not depend on G. However, we can not here determine the time variation of c 
and A completely. 

Now if A = and k = 0, A > one has the constraint 

m + 2n + 3>0. (16) 

Applying the above constraint in eq.(5) one obtains 

R = B'f^d^ j (17) 

where B' = constant. With the aid of eq.(16), eq.(17) implies that 

R > (18) 

as long as 

m + 2n>-2. (19) 



Thus, we see that with a variable speed of light one can solve the horizon and flatness problems of the 
early universe; and justify the present cosmic acceleration. Since our model does not determine m and n 
uniquely, their values has to be found from the present observational data. 



5 Conspirative mechanism (I) 

One can derive the conservation equation from eqs.(l) and (2). This is done by differentiating eq.(l) an 
eliminating R between eqs.(l) and (2). We treat here c and A as variable but G = constant. Consequently, 
one obtains the following equation 

R ( p\ Ac 2 Acc 3kcc , . 

We assume here that the usual energy conservation hold so that 



provided that the other scalars (c, A) conspire to satisfy it, i.e., 

Ac 2 Acc 3/ccc 



8vrG + 4vrG AmGR 2 ° ^ 

In this case one assumes, not only Friedmann equation to be the same for variable c and A, but also the 
usual energy conservation to hold too. With this assumption one can solve eq.(22) to study the effect of 
the variation of c on this cosmology. Again assume the speed of light to have the form 

c = c R n , (23) 

where Co and n are undetermined constants. Integrating eq.(21) using eq.(3) to obtain 

p = A'R- 3 "/ , (24) 

where A' = constant. Integrating eq.(22) using eq.(23) one gets 

A = . 3 " . R- 2 + FR~ 2n , F = constant (25) 
(n- 1) 



If one sets F = 0, we get 



A-^*-«, (26) 
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a variation law that already suggested by several authors (Ozer & Taha, 1986, 1987; Chen & Wu, 1990) 
with different motivation. It is interesting to link the variation of A to the variation of the speed of 
light. We observe that if n = then A = 0. This implies that if the speed of light is constant then 
the cosmological constant vanishes. However, if k = then A = whether c is constant or not. Thus 
one can solve the cosmological constant and flatness problems simultaneously, in radiation and mater 
dominated phases. These problems have preoccupied many scientist for long time. Some string theorists 
believe that A = according to string theory. We therefore, provide a viable model realizing this belief. 
One then would attribute a non-zero cosmological constant to a variation of c. Cosmologists have been 
seeking some symmetry that dictates A to vanish. The existence of a non-zero A in the early universe is 
needed for inflationary solution. The present observations coming from type la supernovae suggest that 
our universe is accelerating at the present time. With the present model we have shown that this cosmic 
expansion is also featured. Thus our model could give a plausible answer to this enigma. However, we 
have shown that the cosmological problems that inflation thought to solve can be solved in the framework 
of the present cosmology. We come to the conclusion that a vanishing cosmological constant follows from 
either the constancy of the speed of light and/or flatness of the universe. 

Again, this setting solves the horizon and flatness problems in the early universe with the same condition 
as the one in eq.(6). Substituting eqs.(24) and (26) in eq.(l) one gets 

This time the flatness and horizon problems can be solved with 

n = -\ or n <^(2-3 7 ). (28) 

For the radiation (7 = |) and matter (7 = 1) dominated phases, we obtain the usual Einstein-de Sitter 
solutions, viz., 

Rcct^ , Rotti , A = 0. (29) 

6 Variable c, A and G 

We will consider here, for completeness,, a model in which c, G and A vary with scale factor (R). 

Accordingly, eqs.(l) and (2) yield (by differentiating eq.(l) and substituting it eq.(2) to eliminate R, R) 

. 3 / p\ G' 3kcc' A'c 2 Acc' , , 

P+ R [ p+ ^) +P G-^ + ^G + ^G=^ (30) 



where the prime ' denotes derivative with respect to the scale factor (R). A simplest solution of eq.(30) 
is to consider the power law for c, A and G of the form 



c = c R n , A = A R m , G = G R a . (31) 
Applying eq.(31) to eq.(30) using eq.(3), one obtains 

P = ART*-* + f ^nfc - nA + A 1 2ra . Q _ 2 

H 4vrG (2n + 37-2) ' v ; 

where for consistency we take m = —2, and A= const. Substituting eq.(32) in eq.(l) using eq.(31) we get 

(R\ 2 = SttGqA 37 c§(2fc-37fc + 7 Ao) p2n-2 f33) 
l^/ 3 (2n + 37-2) ' 1 ' 

The horizon, flatness and cosmological constant problems will be solved with the condition 

n<l(2-3 7 ), (34) 

which is the same as eqs.(6) and (28) (except here m = —2). Equations (32) and (33) generalizes eqs.(8) 
and (9) of Barrow & Magueijo (1999) in their attempt to solve the quasi-flatness and quasi-lambda 
problems. Thus, whether G varies or not, the solution of the horizon and flatness problems is not 
affected. This solution depends only on the variation that A and c assume. This is also evident as a 
does not enter in the constraint in eq.(34), and from eq.(33). The effect of a positive a is to increase the 
energy density of the universe in the early times and decrease (dilute) it at late times. This may help 
explain why, today, we observe the universe to have a deficit in its anticipated energy density (as favored 
by inflationary models). Hence, one need not assume any dark matter to scale up the energy density to 
the required level. A negative a will have the opposite contribution. For a more realistic model of the 
universe, one should consider the effect of bulk viscosity (since it is a basic property of any real fluid) on 
the evolution of the universe. In an earlier work (Arbab, 1997), we have shown that in a viscous universe 
with variable G and A a de-Sitter inflationary solution arises naturally without being imposed. We will 
tackle this problem in the forthcoming letter. 

6.1 Radiation and Matter dominated phases 

We will study here the following two cases: 
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6.1.1 .4 = . 

To solve eq.(33) for matter dominated phase we set 7 = 1, so that for Aq > , one obtains 



R = N m tTT^J , Nm = c J ^L- k A ; n ^\,~. (35) 



This solution is also obtained by Barrow and Magueijo (1999) with G constant [see eq.(15)]. 

4 

3' 



To solve eq.(33) for radiation dominated phase we set 7 = I, so that that for Aq > , one obtains 



R = N r tTT^ , N r = coJ ^° 3*0 ^ n/ i,_i. (36) 

y 3(n + l) 

We see that the scale factor R has the same form but n assumes different values in the two eras (RD & 
MD). This really shows that the way c varies affect the expansion of the universe that result in solving 
may of the cosmological problems. 

6.1.2 A ^ . 

In this case we will assume that 

(2k - 3-fk + 7A0) = , and 2n + 37 - 2 ^ , (37) 

in both eras (RD & MD). This is guaranteed if one takes k = and Ao = 0. For this situation the 
solution of eq.(33) reduces to the Einstein -de Sitter solution, viz., 

Rozt^ (RD) and Recti (MD). (38) 

7 Conspirative mechanism (II) 

We assume again that the usual energy conservation 



to hold. Integrating eq.(39) one obtains 

p = R-^ A = const. (40) 

Using eq.(39) in eq.(30) one gets 

G 1 3kcc' A'c 2 Acc' n 

p g~^gW + ^g + ^g=^ W 



As before we adopt the power law for c and A of the form 



c = c R n , A = A R m . (42) 

Using eqs.(42) and (40), eq.(41) can be integrated to give 

^ Cn (3nk — nAn + An) ~„ 9 „ „ 

G = — -^R 2n - 2 + B , m = -2 , £ = const. (43) 

4vrA (2n + 37-2) ' ' v ; 

Substituting eqs.(40), (42) and (43) in eq.(l) yields 

(R\ 2 cg( 7 A -3fc 7 + 2A;) 87rAi? 

l^J " (2n + 3 7 -2) * ' jD -^ _ - (44) 

We observe that 

A oc R~ 2 , (45) 
a variation that is favored by several cosmologists (e.g., Chen & Wu, Ozer & Taha). 
7.1 Radiation and matter dominated phases 
We study here two cases as in 6.1: 



7.1.1 D=0 . 

In this case eq.(44) reads 

(R V _ C 2( 7 A - 3/c 7 + 2k) R2n _ 2 ^ 



i? / (2n + 3 7 - 2) 



This can be integrated to give 



i 



B-X,„l" ■ N m = coJ ( ( ^° + ^ , n^l,~. (47) 



during matter dominated phase, and 



R = N r W-") , AT r = coy ( 2 3 |° + ^ , 1,-1, (48) 

during radiation dominated phase. We notice that eqs.(35) and (36) are identical with eqs.(46) and (47), 
respectively. This shows that our assumption of eq.(39) was correct. We have already applied a similar 
conservation law for a bulk model (Arbab, 1997). 
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7.1.2 D^O. 

In this case we will assume that 



(2k - 3-fk + 7 A ) = , and 



2n + 37 - 2 ^ 



(49) 



in both eras (RD & MD). This can be satisfied if one takes k = and Ao = 0. For this situation the 
solution of eq.(33) reduces to the Einstein -de Sitter solution, viz., 



Which is a solution we have already obtained in eq.(38). Thus, in an evolving universe particles interact 
(conspire) in such a way the energy conservation law holds. Hence, energy conservation is indeed a 
stringent law that governs the dynamic of the interacting particles in our expanding universe. 

8 Avelino & Martins cosmology 

Avelino & Martins (1999) propose a generalization of general relativity that includes the variation of 
c and G. They however, argued that the this theory is both covariant and Lorentz invariant. This is 
parameterized by writing 



i?oc£2 (RD) and R oc t§ (MD). 



(50) 



= const = A . 

& 

We apply this equation in eq.(20) with the variation 



(51) 



C = CqR 



1 n 



(52) 




R~ 2 + MR- 2n -^ , 



M = const., 



(53) 



Now substitute eqs.(51)- (53) in (1) to get 




D' = 



3 



(54) 



It is clear that this solves the flatness and horizon problems if 



n<-(2-3 7 ) , 



(55) 
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which is the same as eq.(34). However. Avelino and Martins claim that the horizon and flatness problems 
require similar conditions to the one found in the context of the standard cosmological model. They argued 
that a theory that reduces to General Relativity in appropriate limit and solves the horizon and flatness 
problems of the standard model must either violate the strong energy condition or Lorentz invariance, or 
covariance. They also show that the variation in eq.(51) guarantees conservation of mass energy density 
p. The solution of eq.(54) is already worked out in sec. 5, 6, 7. We notice that Barrow (1999, 2003) 
has obtained many solutions that we have found in this letter by different methods. We conclude that 
Avelino & Martins cosmology is equivalent to Albrecht and Magueijo, Barrow, and Moffat models. 

9 Conclusion 

We have studied here the effect of the variation of speed of light and cosmological constant in cosmology. 
We have found that such a variation can in principle solve many of the long persisting problems of 
standard model of cosmology. Of these problems are the flatness, horizon, monopole and cosmological 
constant problems. A conspiracy between these variable constants can still preserve some of the standard 
model picture. The variation of the gravitational constant is not dealt with here, but will be one of our 
future endeavors. These models share some features of the original Moffat, Albrecht and Barrow scenario. 
Our model provides a good solution to the existence of dark matter in the universe if the gravitational 
constant varies with time. It also helps understand the cause of the present cosmic acceleration. Moreover, 
the cosmological constant is found to relax as A oc R~ 2 which is suggested by some cosmologists. 
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